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Abstract. The line graph CG of a directed graph G has a vertex for every 
edge of G and an edge for every path of length 2 in G. In 1967, Knuth used 
the Matrix- Tree Theorem to prove a formula for the number of spanning trees 
of CG, and he asked for a bijective proof [6]. In this paper, we give a bijective 
proof of a generating function identity due to Levine [7] which generalizes 
Knuth's formula. As a result of this proof, we find a bijection between binary 
de Bruijn sequences of degree n and binary sequences of length 2 n ~ 1 . Finally, 
we determine the critical groups of all the Kautz graphs and dc Bruijn graphs, 
generalizing a result of Levine [7|. 



1. Introduction 

In a directed graph G — (V, E), each edge e € E is directed from its source s(e) 
to its target t(e). The directed line graph CG of G with vertex set E, and with an 
edge (e, /) for every pair of edges in G such that t(e) = s(f). A spanning tree of G 
rooted at a vertex r is an edge-induced subgraph of G in which there is a unique 
path from v to r, for all v G V. 

We denote the indegree and outdegree of a vertex v by indeg(w) and outdeg(v), 
respectively, and we denote the number of spanning trees of G by k(G). Knuth 
proved that if every vertex of G has indegree greater than 0, then 

k(CG) = k{G) Y[ outdeg(w) indos(,;) - 1 

v£V 

Knuth's proof relied on the Matrix- Tree Theorem. In his paper, he noted that the 
simple form of this result suggested that a bijective proof was possible, but that it 
was not at all obvious how to find such a bijection [BJ. 

In fact, there are even stronger relations between k(£G) and k(G). Let S 
V} and {x e \e G E} be variables indexed by the vertices and edges of G. The vertex 
and edge generating functions of G are defined as follows, where the sums are taken 
over all rooted spanning trees T of G. 

T eeT T eET 

Levine used linear algebraic methods to prove the following generalization of Knuth's 
result. Our first result in this paper is a bijective proof of Levine's theorem, which 
yields a bijective proof on Knuth's theorem as a special case. 
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Theorem 3.1 Let G — (V, E) be a directed graph in which every vertex has 
indegree greater than 0. Then 

indeg(v) — 1 

n vertex (CG) = n edae (G) II I E Xe 

vEV \s(e)=v 

Using this bijection, we are able to answer the following open question posed by 
Stanley. 

Exercise 5.73 from [8\. Let B(n) be the set of binary de Bruijn sequences of 
degree n, and let S n be the set of all binary sequences of length 2™. Find an ex- 
plicit bijection B(n) x B(n) — ► S(n). 

The critical group K(G) of a graph G is a finite abclian group whose order is 
the number of spanning trees of G. Critical groups have applications in statistical 
physics |4J , algebraic combinatorics |7j , and arithmetic geometry p~] . We review the 
definition of this group in section [2j 

The Kautz graphs Kautz„(m) and the de Bruijn graphs DB n (m) are families of 
iterated line graphs. Kautzi(m) is the complete directed graph on m + 1 vertices, 
without self-loops, and DBi(m) is the complete graph on m vertices, with self- loops. 
These families are defined for n > 1 as follows. 

Kautz„(m) = r^Kautz^m), DB n {m) = C n ^ 1 DB 1 (m) 

Levine recently determined K(DB n (2)) and if(Kautz„(m)), where m is prime [7]. 
We generalize these results, proving the following characterizations of the critical 
groups of all the Kautz and de Bruijn graphs. 



Theorem 5.2 The critical group of DB n (m) is 

K(DB n (m)) = (Z m „) m ~ 2 8 (Z m< ) 



m - 2 (Ts ST\ (7/ \m° _1_ '(m-l) 2 



Theorem 5.3 The critical group of Kautz n (m) is 



n-2 

ll-ifBffl ,\ m2 - 2 en £TS (7/ \ni n - 2 - z (m-l) 2 (m+l) 



X (Kautz„(m)) = (1 m+ i) m ^ © (Z ro »-i) m © (Z^)' 



»=i 



The rest of this paper is organized as follows. In Section |2] we provide background 
and definitions. In Section [3] we introduce a bijection which proves Theorem |3.1| 
We apply this bijection in Section [4] to construct a bijection betweeen binary de 
Bruijn sequences of order n and binary sequences of length 2™ _1 . Finally, in Sec- 
tion [5] we prove Theorem 5.2 and 5.3 giving a complete description of the critical 



groups of the Kautz and de Bruijn graphs. 



2. Background and definitions 



In a directed graph G = (V, E), each edge e £ E is directed from its source s(e) 
to its target t(e). 
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Definition 2.1 (Directed line graph). Let G = (V,E) be a directed graph. The 
directed line graph CG is a directed graph with vertex set E, and with an edge 
(e, /) for every pair of edges e and / of G with t(e) = s(f). 




Figure 2.1. A directed graph and its line graph. 

At times we may speak of a subset F of E as a subgraphs of G - in this case we 
mean the subgraph (V, F). If H is a subgraph of G and v is in H, we denote the 
indegree of v in H by indeg ff (?;), and the outdegree by outdeg ff (w). 

Definition 2.2 (Oriented spanning tree). Let G — (V, E) be a directed graph. An 
oriented spanning tree of G is an acyclic subgraph of G with a distinguished node, 
the rooi, in which there is a unique path from every vertex v £ V to the root. We 
refer to these trees as spanning trees. 

Let T be a spanning tree of G. Every vertex of G has outdegree 1 in T, except 
the root, which has outdegree 0. We denote the number of spanning trees of G by 
k(G), and the number of spanning trees rooted at r by k(G, r). 

Let G = (V, E) be a strongly-connected directed graph, and let Z y be the free 
abelian group generated by vertices of G - the group of of formal linear combinations 
of vertices of G. We define A„ £ 1 V , for all v £ V, as follows. 

a v = J2 (*( e ) - w ) 

e&E s.t s(e)— 1> 

The sandpile group K{G, r) with smfc r is the quotient group 

K(G,r) = Z v /(r,A v \v £ V\r) 

It is well-known that the order of K(G,r) is k(G, r). 

A directed graph G is Eulerian if indeg(w) = outdeg(i;) for all vertices v in V. 
According to Lemma 4.12 of |4 , if G is Eulerian, the sandpile groups if(G, n) and 
JsT(G, r2) are isomorphic for any two 7*1, ^2 in 1/. In this case, we call the group the 
critical group K(G). 

Definition 2.3 (The Laplacian). Let G = (V^E 1 ) be a finite directed graph with 
vertices V\,V2, ■ ■ .viv'l- The adjacency matrix A(G) of G is the \V\ x \V\ matrix in 
which A(G)ij is the multiplicity of the edge (Vi,Vj) in E. The degree matrix D(G) 
is the I V| x \V\ diagonal matrix in which Da = outdeg(fi). The Laplacian L(G) of 
G is defined as A(G) — D(G). 

Note that the row vectors of L(G) are the elements A„. We consider L(G) T as 
a Z-linear operator on 1 V - its image is the subgroup generated by the A„. For a 
strongly-connected Eulerian graph G, the Laplacian has exactly one eigenvalue 0, 
so for such a graph G, we have 

Z^/im L{G) Z y /im L(G) T Z © K(G) 

The following elementary row and column operations on matrices with entries in a 
ring i? are invertible over i?. 
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- Permuting two rows (columns) 

- Adding a multiple of a row (column) by an element of R to another row (column) 

- Multiplying the entries of a row (column) by a unit 

If V is obtained from L(G) by invertible row and column operations over Z, then 
Z y /Im V = Z v /Im L(G). 

Suppose that R is a principal ideal domain. Under these operations, any matrix 
with entries in R is equivalent to a matrix in Smith normal form. A matrix in this 
form is diagonal, and its diagonal entries in, X22, ■ • ■ x nn are elements of R such that 
x (i+i)(i+i) is a multiple of xu for all i < n. These entries are called the invariant 
factors of the original integer matrix, and they are unique up to multiplication by 
units. If the invariant factors of L(G) over Z are xu, ■ ■ ■ x nn then 

n 

Z v /lm L(G) = 0Zx„ 
t=i 

Thus, row-reducing the Laplacian yields information about the critical group. 



3. Counting spanning trees 

Let G = {V,E) be a directed graph, and let {x v } v ^v an d {Xe}eeE be variables 
indexed by the vertices and edges of G. The edge and vertex generating functions, 
which enumerate the spanning trees of G, are defined as follows 



n^(G) = J2U 



T e<ET 



K V£rteX (G) = V TT 



^ 11 ~t(e) 
T eST 

where T ranges over all spanning trees of G. In this section, we give a bijective 
proof of the following identity, solving a problem posed by Levine in [7 

Theorem 3.1. Let G=(V,E) be a directed graph in which every vertex has indegree 
greater than 0. Then: 

(\ indeg(v) — l 
s(e)=v I 

In order to find a bijection, we adopt the following strategy. We put an arbitrary 
total order on the edges in E. 

- We provide a bijection between monomial terms on the right-hand side of Eq. 



(3.1 1 and tree arrays, which are arrays of lists, one list for each vertex v G V. 



Then we present a map a that take a tree array to a spanning tree of CG which 



contributes the same term to the left-hand side of Eq. (3.1 1. 
- Finally, we show that a is bijective by constructing an inverse map tt which takes 
a spanning tree of CG to a tree array. 

We define a list to be an ordered tuple of edges. We append an element x to a 
list I by adding x to the end of We pop list I by removing the first element of I. 
We denote the number of times an clement e appears in a list I by N(l, e) . 

Let v be a vertex of G and let l' v be a list with indeg(v) — 1 elements, all of which 
are edges with source v. We map l' v to a monomial term of (J2 s (e)=v x e ) lndce ^'~ 1 , 
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as follows. 

l v e 2> ■ • ■ Gindcg(v) — l) y %ei%t 

This map provides a bijection between lists l' v and terms of Q2 s (e)=v £e) mdog ^ t '' l_:L - 



Therefore, a term on the right-hand side of Eq. (3.1l corresponds to a choice of 



spanning tree T of G and a choice of one such list l' v for each vertex v. 



Suppose a monomial term on the right-hand side of Eq. (3.1 ) corresponds to a 
spanning tree T rooted at r and an array of lists (l' v ). For each vertex v £ V\r, 
we obtain l v by appending the unique edge e in T with source v to the list V v . We 
obtain l r by appending a new variable Q, to V r . 

Each list l v has length indeg(i>), for v £ V. We call an array of lists (l v )vev 
obtained in this way a tree array. By construction, terms on the right-hand side of 



Eq. (3.1) are in bijection with tree arrays. 

We now define the bijective map a, which takes a tree array of G to a spanning 
tree of CG. 

The bijection a: We start with a tree array (l v ) and an empty subgraph T" of 
CG. Then we run the following algorithm. 

Step 1. Let R be the subset of edges e of G for which N(l s r e -\, e) = and outdeg T /(e) = 
0. Let / be the smallest edge in R under the order on E. 

Step 2. Pop the first element g from the list lt(f)- If <? is ^2, then a((l v }) = J". 

Step 3. Otherwise, g £ E and s(</) = £(/)■ Add the edge (/, g) to T', and then 
return to step 1. 

We also define a map ir which takes a spanning tree of CG to a tree array of G. 

The inverse map ir: We start with a spanning tree T" of CG, and an empty list 
l v at each vertex v £ V. This map is given by another algorithm. 

Step 1. Let S be the set of leaves of T'. Let / be the smallest edge in S under the 
order on E. 

Step 2. If / is not the root of T', remove / and its outedge (f,g) from T' , and 

append g to h(f)- Go back to step 2. 
Step 3. If / is the root of I", append £1 to hlf), and return the array of lists. 

As an example, we apply a to a tree array in a small directed graph G. We order 
the edges of G by the lexigraphic order. 

li=[0] 
l 2 =[(2,5)] 
l 3 =[(3,4)] 
l 4 =[(4,l),(4,l)] 

X25X34X41X54X41 l 5=[( 5 - 4 )] 

Figure 3.1. The graph G, with a spanning tree T highlighted in red. Below the 
graph is a monomial term of n vertex (G) , where Xij is the variable for edge (i,j)- 
The tree array corresponding to this term is shown to the right. In the term and 
the tree array, red elements correspond to edges of the tree. 
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© ll=[Q] © ll=[0] 

\ 4 ®\ W l 3 =[(3,4)] \t®\^ h=l&rf] 

®s } l 4 =[(4,l),(4,l)] ® St y | 4 =[(4,1),(4,1)] 

© l 5 =[(5,4)] © l 5 =[(5,4)] 

Figure 3.2. The first two edges added to T' by the algorithm for a. Initially, the 
edges (1,2) and (1,3) do not appear in the lists. We pop (2,5) from I2 and add the 

edge ((1,2), (2,5)) to T' . Then the edges (1,3) and (5,4) have outdegree in T' 
and do not appear in the lists. We pop (3,4) from I3 and add ((1, 3), (3,4)) to T' . 



© 



© 



r® 

J' 

! ® 



li=[Q] 

l 4 =[(4,l),(4,l)] 



© 



* i © 



© 



li=[0] 

i 2 =[^5T] 
l 3 =[pj<] 
l4=U*rtf,(4,l)] 




li=[Q] 
l 2 =[Pr5T] 
l 3 =[pr<] 

i 4 =[j^T,^rT] 

l 5 =U5^] 

T/ie Zasi i/iree edges added to T' , and the final tree, 
left in the lists of the tree array is £1. 

In order to prove Theorem |3.1[ we first prove three lemmas. In the definition of 
the algorithm for the map a, we assumed that the set R is always non-empty in 
step 1 and that the the list U(f) is always non-empty in step 2. In Lemma 3.2 we 
show that both assumptions are valid. 



li=[0l 

l 2 =[Pr5T] 
l 4 =[^rfT,^rlT] 

T/ie /as£ element 




Lemma 3.2. TTie algorithm used to define map a is well-defined: at step 1, the set 
R is non-empty, and at step 2, the list lt(j) is non-empty. 

Proof. After k edges have been added to T', there are \E\ — k elements left in all 
the lists l v , where one of the elements is il. There are \E\ — k — 1 edges left in the 
lists, but there are \E\ — k edges of G which do not have an outedge in T' , so R 
must be non-empty in step 1. 

Every time we pop l v , we add an edge (/, g) to T", where t(f) = v. When we are 
at step 2, outdeg T ,(/) = 0, so at most indeg(i(/)) — 1 of the elements of l t (f\ have 
been popped. Therefore, the list lt{f) is always nonempty at step 2. The algorithm 
is well-defined. □ 



The following lemma shows that a takes a tree array corresponding to a term 



on the right-hand side of Eq. (3.1 1 to a spanning tree which contributes the same 
term to the left-hand side. 

Lemma 3.3. Suppose that (l v ) is a tree array and that o~((l v )) = T' . Then T' is a 
spanning tree of CG, and indeg T ,(e) — N(l s ^,e), for all ee£. 

Proof. Let 1(e) be the initial value of N(l s ( e ),e). By the definition of a tree array, 
the edges which are the last elements of the lists l v form a spanning tree T of G. 
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We claim that T' is acylic, because the last edge of a cycle is never included in 
T' . While the algorithm is running, suppose that (e„, ei) is not an edge of T", and 
that it completes a cycle (ei, e2), (e2, 63), . . . (e„_i, e n ) of edges in T'. Since (ei, e2) 
was already added to T', JV(Z s ( ei ),ei) must be 0. Therefore, (e„,ei) will never be 
added to T . 

We say a vertex we Vis cleared if all the elements of its list are popped. Suppose 
that e = (v, w) is an edge in T. The list is cleared when all the edges of G with 
target to have an outedge in T 1 . Then w can only be cleared after an outedge (e, /) 
of e is added to T. The edge (e, /) can only be added to T" when N(l v ,e) = 0. 
Because e is an edge of T, it is the last element of l v , so v must be cleared before 
to can be cleared. 

The algorithm terminates when ft is popped from l r , which occurs when r is 
cleared. There is a path (t>, t>i, V2, ■ ■ ■ Vk, r) in T from any vertex v to r. Therefore 
r can only be cleared after all the vertices on this path are cleared. Thus, all the 
vertices of G are cleared when the algorithm finishes, so there are \E\ — 1 edges in 
the subgraph T' . 

All the vertices of CG has an outedge in T", except one. Since T" is acyclic, it 
is a spanning tree of CG. Because indeg T , (e) + N(l s ^,e) is constant, when the 
algorithm returns T", indeg T /(e) = 1(e) for all e E E. □ 

In our final lemma, we show that ir will take a spanning tree T" of CG and 
reconstruct a tree array (l v ). 

Lemma 3.4. Suppose T' is a spanning tree of CG with root r' , and that n takes 
T' to the array of lists (l v ). Then (l v ) is a tree array, which means that 

(a) The length of l v is indeg(v), for all v E V. 

(b) Every element of l v is an edge with source v, for all vertices v except t(r'). 
The last element of ' lt( r ') is Q> an d every other element ofl t ( r ') is an edge with 
source t(r'). 

(c) The set T of edges which are the last elements of the lists {l v \v € V\t(r')} is a 
spanning tree of G. 

Proof. We first show parts (a) and (b). Each time an edge e G E is removed from 
T' , an element is appended to the list l t { e )- Since r' can only be removed after all 
the other vertices of T", this algorithm adds indeg(w) elements to l v for all v G V, 
so part (a) holds. Every element of the list l v is an edge with source v, with the 
exception of fi, which is the last element of l r , so part (b) holds. 

While the algorithm tt is running, say a vertex v G V is filled if l v has indeg(w) 
elements. Every vertex is eventually filled, so the order in which vertices arc filled 
is a total order on V. 

We claim that this order is a topological sort of the subgraph T. Suppose that 
/ = (v, to) is the last element of l v for some vertex t; other than r. Vertex v was 
filled at step 2, right after some leaf e and some edge (e, /) were removed from T" 
in step 1. However, to cannot be filled until / is removed from T', which happens 
after e and (e, /) are removed from T'. Therefore v is filled before to, so the filling 
ordering on V is a topological sort of T, and T is acyclic. 

After tt terminates, t(r') has no outedge in T and every other vertex of G has 
one outedge. Then T is a spanning tree of G, and part (c) holds. □ 



We now prove the main result. 
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Proof of Theorem \3.1\ Let (l v ) be a tree array and let T' = a((l v }). We show 
the following claim by induction on n: after n edges have been added to T' by the 
algorithm for a((l v )), and n edges have been removed from T' by the algorithm for 
7r(T'), we have 

(a) The set of edges added to T' by a is the set of edges removed by tt. 

(b) The elements popped from l v by a are exactly the elements added to l v by n, 
in the same order. 

In the base case n = 0, both claims hold trivially. Suppose both results hold for 
n = k. The edge e is a leaf of T' in tt if and only if it satisfies N(l s r e -\, e) = and 
outdeg T /(e) = in a. 

Therefore, the (k + l)st edge (/, g) added to T' by a is also the (k + l)st edge 
removed from T' by tt, and the element g popped from lt(f) in f is also the element 
appended to h(f) by tt. Both claims hold for n — k + 1. By induction, they hold 
for all n<\E\ - 1. 

When n = — 1, condition (b) implies that 7r(T") = Then 7r is a left 

inverse of a, and a is injective. 

By similar reasoning, tt is a right inverse of a, and er is surjective. So cr is 
a bijection between tree arrays in G and spanning trees of CG. The bijection <r 



induces between equal terms in Eq. (3.1l proves Theorem 3.1 □ 



4. The de Bruijn bijection 

A binary de Bruijn sequence of degree n is a cyclic binary sequence B such that 
every binary sequence of length n appears as a subsequence of consecutive elements 
of B exactly once. For example, 0011 is a binary de Bruijn sequence of degree 2, 
since its cyclic subsequences of length 2 are 00, 01, 11, and 10. 

It is well-known that there are 2 2 binary de Bruijn sequences of degree n. 
Stanley posed the following open problem in |5]. 

Exercise 5.73 of [5]. Let B(n) be the set of binary de Bruijn sequences of degree 
n, and let S n be the set of all binary sequences of length 2™. Find an explicit 
bijection B(n) x B(n) — * S n . 

Our solution to this problem involves the de Bruijn graphs, which are closely 
related to de Bruijn sequences. 

Definition 4.1 (de Bruijn graph). The de Bruijn graph DB n (m) has m" vertices, 
which are identified with the strings of length n on m symbols. The edges of the 
graph are labeled with the strings of length n+1 on 7H symbols. The edge soSi • ■ ■ 
has source sqSi . . . s n _i and target S1S2 ■ ■ ■ s n . 

An edge of DB n (m) can be identified with the vertex of DB n+ i(m) that is 
labeled with the same string of length n+1. With this identification, we have 

DB„ (to) = CDB n -i{m) 

Each vertex v = soS\...s n -i of DB n {2) has two outedges, s\S2, ■ ■ ■ s n _i0 and 
S1S2 • • • s n _il. We call these edges the zero edge of v and the one edge of v, 
respectively. 

It is well-known that binary de Bruijn sequences of degree n are in bijection 
with Hamiltonian paths in DB n (2). Let B = b bi . . . &2 n -i be a binary de Bruijn 
sequence of degree n. Let «j = bib i+ i . . . &; + „_i, for < i < 2 n — 1, where indices 
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are taken mod 2" . The path (vo, V\, . . . «2"-i) is the corresponding Hamiltonian 
path in DB n {2). 

Theorem 4.2. There is an explicit bijection between B{n) and the set of binary 
sequences of length 2 n ~ 1 , for n > 1. 

Proof. We describe a bijection between Hamiltonian paths in DB n {2) and binary 
sequences of length 2 n_1 . By composing this bijection with the map between de 
Bruijn sequences and Hamiltonian paths, we construct the desired bijection. 

We order the vertices in DBk(2) by the lexicographic order on their associated 
binary strings, for 1 < k < n. Let [v\, . . .v 2 ™) be a Hamiltonian path in DB n (2). 
This path is an oriented spanning tree of DB n (2), so we can apply the inverse map 
7r defined in Section [3] to it. 

Let A n _i be the tree array ^4„_i = 7r(t>i, . . . v 2 ^)- We recursively define a se- 
quence of tree arrays Ak, for 1 < k < n — 1. Suppose we have a tree array A^+i 
in DBk+\{2) . Let Tk+\ be the spanning tree consisting of the edges which are the 
last elements of the lists in Ak+i- We define Ak to be n(Tk+i). 

We construct a binary sequence sis 2 ■ ■ ■ s 2 n-i from these tree arrays. We denote 
vertex w's list in the tree array Ak by (Ak) w - Let S2™- 1 be if the first element of 
{A n _i) s r V2n \ is the zero edge of s(v 2 ™), and 1 otherwise. 

We define s 2 k through s 2 *>+i_i, for 1 < k < n — 2, as follows. Let u>i, w 2 , ■ ■ ■ w 2 k 
be the vertices of DBk(2), in lexicographic order. Let S2 fc +i-i be if the first 
element of (Ak) Wi is the zero edge of w i: and 1 otherwise. Let si be if T\ is rooted 
at vertex 0, and 1 otherwise. 

The string S\S-2 ■ ■ ■ s 2 n-i is the binary sequence that corresponds to the Hamil- 
tonian path we began with. 

Now we construct the inverse map, from binary sequences to Hamiltonian paths. 
Given any binary sequence S of length 2™ _1 , we use the first 2™ _1 — 1 characters of 
the sequence to invert the previous procedure and construct a sequence of spanning 
trees Ti,T 2 , . . .T n _\. The tree Tj, will be a spanning tree of DBk(2). 

We determine Tk recursively. The tree T\ in DBi(2) is rooted at if si is 0, 
and rooted at 1 otherwise. Assume that the first 2 fe — 1 characters of S determine 
a spanning tree Tk of DBk(2), where k < n — 2. We choose a tree array Ak of 
DBk(2) using this tree and the next 2 k characters of S, as follows. 

Let the vertices of DBk(2) be w±,W2,-- -w 2 k, in lexicographic order. The first 
element of (Ak) Wi is the zero edge of Wi if S2 fc +i-i is 0, and the one edge of w.j 
otherwise. The second element of (Ak) Wi comes from Tk- We define Tfc+i to be 
a(Ak), using the map defined in Section [3j 

We use T n —i to construct a tree array A n -\ such that a{A n -i) is a Hamiltonian 
path in DB n (2). Let r be the root of T n _i, and let v be another arbitrary vertex. 
The list l v in the array A„_i must contain two distinct edges, if a(A n _i) is a 
Hamiltonian path. The second edge in l v must be the unique edge in T n _i with 
source v, so l v is determined. Our only remaining choice is which of the two edges 
of DB n -\(2) with source r to include in l r , which we determine by s 2 n-i. 

Clearly, this map from binary sequences to Hamiltonian paths inverts the map 
from Hamiltonian paths to binary sequences. Therefore, our first map is the bijec- 
tion we need. □ 
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This bijection can easily be generalized to count the k-ary de Bruijn sequences, 
in which the 2-symbol alphabet {0,1} is replaced with the fc-symbol alphabet 
{0,1,... k -I}. 

5. The Kautz and de Bruijn graphs 

In this section, we determine the critical groups of all the Kautz graphs and the 
de Bruijn graphs. The critical groups of these graphs have been found in some 
special cases by Levine [7]. 

The Kautz graphs are similar to the de Bruijn graphs, except that the vertices 
are indexed by Kautz strings. A Kautz string is a string in which no two adjacent 
characters are the same. 

Definition 5.1 (Kautz graph). The Kautz graph Kautz„(m) has (m + l)m™ -1 
vertices, identified with the Kautz strings of length nonmll symbols. The edges 
of the graph are labeled with the Kautz strings of length n + 1 on m + 1 symbols, 
such that the edge sqSi . . . s n has source so s i ■ ■ • s n -i and target s\S2 ■ ■ ■ s n . 

We also consider the Kautz and de Bruijn graphs as families of iterated line 
graphs. Kautzi(m) is the complete directed graph on m + 1 vertices, without self- 
loops, and DB\(m) is the complete directed graph on m vertices, with self- loops. 
Then for n > 1, we have 

Kautz„ +1 (m) = £Kautz„(m) = /T'Kautz^m) 

DB n+l (m) = CDB n (m) = C n DBx{m) 

We say a directed graph G = (V, E) is balanced k-regular if indeg(u) = outdeg(w) = 
k for all v G V. Both Kautz„(m) and DB n (m) are balanced m-regular, for all 
nffN, which implies that they are Eulerian. Since these graphs are also strongly- 
connected, their critical groups are defined. 

Levine found the critical groups of the de Bruijn graphs DB n {2) and the Kautz 
graphs K n (j>), where p is prime [7 . In this section we characterize the critical 
groups of all the Kautz and de Bruijn graphs. We prove the following theorems. 



Theorem 5.2. The critical group of DB n (m) is 

n-l 

K (DB n (m)) = (Z m ~) m ~ 2 (Z r 

i=l 

Theorem 5.3. The critical group of Kautz n (m) is 

K(^ Zn (m))^(z m+1 r- i ©(z m „- 1 r 2 " 2 ®0(z m .r' J "- i|j "" +u 

i=l 



(m-1) 2 



In order to prove these theorems, we first prove two lemmas about row-reducing 
the Laplacians L(Kautz n (m)) and L(DB n (m)). We refer to the row and column 
of a vertex v in the Laplacian by R(v) and C(v), respectively. We also use L(v,w) 
to denote the entry in the row of v and the column of w. 

We say two strings of length n are similar if their last n — 1 characters are equal. 
Similarity is an equivalence relation. We partition the vertices of Kautz„(m) and 
DB n (m) into equivalence classes, by grouping vertices labeled with similar strings 
in the same class. There are m vertices in each class. 



COUNTING THE SPANNING TREES OF A DIRECTED LINE GRAPH 



11 



Lemma 5.4. Let G = (V, E) be a Kautz graph Kautz n+ i{m) or a de Bruijn graph 
DB n+ i(m), where n G N. Then G contains a cycle (v\, v%, . . . v c ) of length c — 
\V\/m which contains one vertex from each class. 

Proof. Let G' , the predecessor of G, be Kautz„(m) if G is Kautz rl+1 (m), and 
DB n (m) if G is DB n+1 (m). 

First we show that there is a Hamiltonian cycle in G' . Such a cycle exists in the 
complete graphs K m and K m+ i, so the case n — 1 is done. There is an Eulerian 
tour of Kautz„_i(r7i) and of D-B„_i(m) for n > 1, since graphs in both families 
are Eulerian. Because G' is either £Kautz„_ 1 (m) or CDB n _i(m), one of these 
Eulerian tours induces a Hamiltonian cycle in G", for n > 1. 

The Hamiltonian cycle in G" can be represented as a string S = s\S2 ■ ■ ■ s„+ c _i, 
where the ith vertex of the cycle is labeled with s^Si+i . . . s;+„-i. 

We use string S to find a cycle in G. Let vi — SiS.; + i . . . Si +n for i < c, and let 
v c = s c ,s c+ i . . . s„4- c _iSi. By the construction of S, (vi, i>2, • ■ ■ v c ) is a cycle which 
contains one vertex from each class. □ 



In the next lemma we show that every invariant factor of L(Kautz„ +1 (m)) and 
L{DB n (m)) is either a multiple of m or relatively prime to m. We prove this lemma 
by row-reducing the Laplacian in an order derived from the cycle in Lcmma|5.4| 



Lemma 5.5. Let G = (V, E) be a Kautz graph Kautz n+ i(m) or a de Bruijn graph 
DB n+ \(m), where n G N. The first c = \V\/m invariant factors of L(G) are 
relatively prime to m, and all of the rest are divisible by m. 

Proof. We reduce the Laplacian L(G) over the principal ideal domain Z m . Let the 
invariant factors of L(G) over Z be x\, X2, ■ ■ ■ x\y\. Any invertible row or column 
operation over Z descends to an invertible operation over Z m , so the invariant 
factors of L{G) over Z m are the x^ mod m. 

Let (vi,V2, ■ ■ ■ v c ) be the cycle in G from Lemma 5.4 and let be the set of 



vertices in the class of i>;. We take indices mod c, so v c+ i is v\. 

Note that if u and v are vertices in the same class, then (u, w) is an edge if and 
only if (u, w) is, for all w. Therefore, the rows of u and v in the adjacency matrix 
A(G) are the same. 

Because every vertex of G has outdegree m, L(G) = A(G) mod m. Therefore 
rows R(u) and R(v) are congruent mod m if u and v are in the same class. We 
reduce the Laplacian in c stages. In the ith stage, we subtract row R(vi) from R(v) 
for all v G After this operation, every entry of R(v) is divisible by m. 
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Figure 5.1. Reducing L(Kautzi(^)) using the cycle (01,12,20). The original 
Laplacian is on the left. We obtain the reduced Laplacian on the right by 
subtracting R(0l) from i?(21), i?(12) from R(02), and R(20) from R(10). Every 
entry of rows i?(02), i?(10), and i?(21) is divisible by 2. 
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The entry L(vi, Wj+i) is 1 before and after these row operations, for 1 < i < c. 
We claim that in the reduced Laplacian, every entry of C{vi + i) is divisible by to 
except L(vi, t>i+i). There are to edges with target Vi + i in G. The sources of these 
edges are the to vertices in [m]. After the row operations, every entry of R(v) is 
divisible by to for v S so L(vi, Uj+i) is the only entry in C{vi+\) which is 

non-zero mod to, for 1 < i < c. 

By permuting rows and columns, we move L(w i ,w i+1 ) to the zth diagonal entry 
La of the Laplacian. The reduced Laplacian is now in the form 







mod 



where I c is the c x c identity matrix. Using column operations, we can make all the 
entries in A divisible by m, without changing the rest of the matrix mod to. After 
this column operations, the Laplacian is in Smith normal form. 

The first c invariant factors of L(G) over Z m are 1, so the first c invariant factors 
of L(G) over Z are relatively prime to to. The last \V\ — c invariant factors of L(G) 
over Z m are 0, so the last |V| — c invariant factors of L(G) over Z are divisible by 
to. The lemma holds. □ 



We use Lemmas |5.4| and |5.5| to characterize the critical group of the Kautz and 
de Bruijn graphs. The first step is finding the orders of these groups. We apply 



Theorem 3.1 to DB n (m), and we let all the variables x e equal 1, to find that 

K(DB n+1 (m)) — k (DB n (m)) (m^" 1 )™") 

The number of spanning trees of the complete graph DBi(m) is to™" 1 4 j. By a 
simple induction, we have 

K(DB n (m)) =m m ™- 1 

In an Eulcrian graph, the sandpilc groups K (G, v) are isomorphic for all vertices 
v, so \K(G)\ = k(G)/|V|. Therefore, we have 

\K(DB n (m)) | = m m "~ n ~ 1 (5.1) 

Similarly, we have 

«;(Kautz n (m)) = (m + i)™ m (™"-D(™+D 

\K (Kautz„(m)) | = (m + l)m-l m (m»+m»- 1 -m-n) (5 2 ) 

We are ready to prove theorems |5.2| and |5.3[ 

Proof of Theorem \5.%\ We proceed by induction on n. The critical group of the 
complete graph on m vertices is (Z m ) m_2 , so the base case holds. 

Assume that Theorem |5.2| holds for n — 1, where n > 1. We prove it for n. As 
shown by Levine [7], if G is a balanced fc-regular graph, then 

kK(CG) Si K(G) (5.3) 

We will use this fact to determine Syl p (K(DB n (m))), the Sylow-p subgroup of 
K(DB n (m)), for any prime p. We break into two cases: either p does not divide 
to, or p divides to. 



If p does not divide to, then by Eq. (5.3 1, we have 



Syl p (K (DB n (m))) = Syl p (K (DB n _i(m))) 
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By the inductive hypothesis, Sylow p (K(DB n (m))) is the trivial group. 

Now let p be a prime that divides to, and suppose p k is the largest power of p 
that divides to. Let the Sylow-p subgroup of DB n (m) be 

SyL, (K (DB n (m))) = T£ © © . . . © Z£j 

By Lemma [5. 5 1 K(DB n (m)) can be written as a direct sum of cyclic groups, such 
that the order of each group is either non-zero mod p or divisible by p k . Thus, 
a, = for i < k. Further, we can derive the order of Syl p (K(DB n (m))) from Eq. 
( phi) . We find that 

I 

Y^iai = k{m n -n-l) (5.4) 

i = n 

because the express ion o n the right-hand side equals the number of factors of p in 
to™ ~ n_1 . By Eq. (5.3 1 and the inductive hypothesis, we know that 

/Syl p (K (DB n (m))) = Z^ 1 © Z^+ 2 © ... © Z£{_» = Syl p (K (DB n _i(m))) 

n-2 

z^+> © z^ 2 © ... © z°u = (z p „ fc ) m ~ 2 © (v) m ""~' (m ~ 1)2 ( 5 - 5 ) 



i=i 

Eq. (5.5) implies that a nk = to — 2, that au + x\k — to"~ 2 ~ 4 (to— l) 2 for 1 < i < n — 2, 
and that <Xj = for i > nfc or k \ i. The only a, which we have not yet determined 
is dfc. We solve Eq. (5.4 1 for o/t, by moving all the other to the right-hand side 
and dividing by k. 

n-l 

a k = (to" -n-l)-n(m-2)-^j (m"" 1- ^ - l) 2 ) 

i=2 

By evaluating the geometric series, we find that a k — m n ~ 2 (m — l) 2 . With these 
values, we may write 

n-l 

Syl p (K(DB n (m)))=(Z pnk ) m - 2 ®($(Z pln ) m "~ 

The Sylow-p subgroups of K(DB n (m)) are trivial for p\m. Taking the direct sum 
of the Sylow-p subgroups over primes p which divide to, we find 

n-l 

K(DB n (m)) S 0Sylow p (K (DB n (m))) S (Z OT „) ro " 2 © (Z mI ) m ""~ l(m ~ 1)2 

p|m i=l 

With this equation we complete the inductive step, as desired. □ 



Proof of Theorem \5.3\ This proof is similar to the proof of Theorcm |5.2| Again, 
we induct on n. Because the critical group of the complete graph on to + 1 vertices 
is (Z TO+ i) m_1 , the base case holds. 

Assume that Theorem 5.3 holds for n — 1, where n > 1. Using Eq. (5.3 1, we 



calculate the direct sum of the Sylow-p subgroups of -R"(Kautz n (m)) over primes p 
which do not divide to, as follows 



Syl p (K (Kautz n (to) ) ) = Syl p (K (Kaut Zl (m))) = (Z m+1 ) r 



(5.6) 
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Now let p be a prime that divides m. Suppose that p k is the largest power of p 
that divides m, and that the Sylow-p subgroup of K (Kautz„(m)) is 

Syl p (K (Kautz„(m))) = 1^ ffi Z a J © . . . © 



Lemma 5.5 implies that ai — for i < k. Furthermore, we know the order of the 
Sylow-p subgroup of K (Kautz„(m)) from Eq. (5.2 1, which implies that 



E 

i—k 



iai = k ( 



m" + m n 1 — m — n) 



(5.7) 



because the expression on t he ri ght-hand side equals the number of factors of p in 



p te Syl p (if (Kautz„(m))) = Z£ fe+1 © Z^+ 2 «... © S Syl p (if (Kautz„_ 1 (m))) 

By the inductive hypothesis, we find that a^+i)k — m n ~ 3 ~' l (m~ l) 2 (m + 1) for 
1 < i < n — 3, that Q( n _i)fc = ™ 2 — 2, and that a.j = for i > (n — l)fc or fc { : 



We solve Eq. (|5.7| for at- We find that au = m 2 — 2 if n 



2 and that 



afe = 77?," — m™ 1 — 777™ 2 + 1 if 77 > 2. Then we have 

2 "~ 2 

Syl p (K (Kautz„(m))) = (Z p( „- 1)fc ) m ~ 2 © (Z pi ,; 



m™ 2 i (m-l) 2 (m+l) 



(5.8) 



By taking the direct sum of Eq. (5.6) and Eq. ( |5.8| over all primes which divide 
m, we complete the inductive step and prove Theorem |5.3| as desired. □ 
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